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Abstract 

We investigate moduh field dynamics in supergravity/M-theory like set ups 
where we turn on fluxes along some or all of the extra dimensions. As has 
been argued in the context of string theory, we observe that the fluxes tend to 
stabilize the squashing (or shape) modes. Generically we find that at late times 
the shape is frozen while the radion evolves as a quintessence field. At earlier 
times we have a phase of radiation domination where both the volume and the 
shape moduli are slowly evolving. However, depending on the initial conditions 
and the parameters of the theory, like the value of the fluxes, curvature of the 
internal manifold and so on, the dynamics of the internal manifold can be richer 
with interesting cosmological consequences, including inflation. 
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INTRODUCTION 



Unified Gravitational theories, dating back to Kaluza-Klein theory till the most mod- 
ern String/M-theory, all seem to require the existence of extra dimensions. Under- 
standing their nature is therefore of key interest. Is the "internal manifold" consti- 
tuted by the extra dimensions stable, i.e. frozen in time, or is it evolving slowly? What 
are the problems and advantages associated with either of these scenarios? Even if 
the internal manifold is stable now, could it have had an earlier dynamic phase? In 
this paper we try to address some of these issues along with the cosmological and 
particle physics implications of these different scenarios. 

Recently in the context of String theory it has been argued that by turning on 
the fluxes one can fix the complex moduli [Q, which loosely speaking corresponds 
to the shape(s) of the internal manifold. One of the motivations for this paper is 
to study this mechanism of stabilizing the shape moduli explicitly in the context of 
supergravity (SUGRA)/M-theory type compactifications |2J using fluxes along some 
or all of the extra dimensions. Hence we start with a general bosonic^ SUGRA type 
action containing the usual Einstein scalar curvature term and various n-form fluxes. 
To keep the analysis as general as possible we also include a higher dimensional 
cosmological constant. For moduli fields, we choose to study the size and the shape 
or equivalently the volume and the squashing mode (preserving a subgroup, say H, 
of the total isometry group G of the internal manifold) which has been previously 
shown to have a consistent dimensional reduction [21 E] . Such consistency arguments 
are known to hold not only for pure gravitional action but also when fluxes are turned 
on along all of the extra dimensions. Here we first show that it is also possible to 
"consistently" turn on fluxes only along a submanifold of the full internal manifold 
which preserves the subgroup H of isometrics. We find, perhaps not surprisingly, that 
these extra fluxes have the effect of stabilizing the shape or the squashing mode in 
congruence with pp. 

Although in general the shape is stabilized, the radion or the volume which is 
analogous to the Kahler moduli is not^; it seems to have a runaway exponential 
potential. The advantage of having such a potential is that the radion can now act 
as a quintessence field jHj, with the potential energy becoming vanishingly small as 
the volume rolls along the exponential. Such a potential can then naturally explain 
the exponential hierarchy between the currently observed cosmological energy density 
< A > and the (reduced) Planck mass Mp jTj 



On the other hand if the radion is stabilized at a potential minima, then one has 
to typically invoke fine tuning to justify why the potential energy vanishes (or is 
very small (jO.ll) ) at the minima^. The flip side of having a runaway radion is that it 
typically couples to matter-radiation. This first of all makes several coupling constants 

^To study the moduli dynamics the fermionic fields can be set to zero. 

^For attempts in this regard see This problem of stabilizing the volume or radion is also 
closely related to the issue of dilaton stabilazation. 

^As to why the quantum corrections to the vacuum energy are also as small as H().l|l is a much 
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like the fine structure and Newton's constant time dependent (as they depend on the 
volume which is now evolving in time) and there are strong constraints jHl CHI on 
the same. Second, since the quintessence field is effectively massless, it mediates a 
fifth force between matter/radiation leading to violation of the equivalence principle 
which has been tested to high precision, leading to strong bounds [TU] . Is it possible 
to avoid these bounds? 

The effective potential that we obtain for the moduli turns out to be a sum of six 
exponentials with complicated but rich features. In particular what is really interest- 
ing is that we find, for certain values of group theoretic parameters and fiuxes there 
are late time attractor solutions where a linear combination of the radion and the 
shape, rather than just the shape, is stabilized! Now note that radiation/matter will 
typically couple to a linear combination of the radion and the shape field. Thus if the 
combination that is stabilized coincides with (or is very close to) the linear combina- 
tion that couples to matter/radiation then the fifth force constraints will be easier to 
satisfy. In our case we find that although the introduction of the extra fiux allows us 
to stabilize a linear combination which is close to the "coupling" linear combination, 
it is not close enough to avoid all the bounds. However, the results are encouraging 
and suggest that it may be possible to obtain the right linear combination by turning 
on fiuxes in other directions or by accounting for potentials coming from brane-gas 
jTT] etc. Also, recently discovered chameleon effects of the scalar moduh fields |T2] 
may play a role in ameliorating these constraints. Both the linear-combination and 
chameleon effects also have a tendency to moderate the constraints coming from time 
variation of physical constants. These latter type of constraints however can also be 
avoided in another way. Since the kinetic energy of the quintessence field is related 
to the potential energy or the cosmological constant through the equation of state 
parameter, one can check that (jO.ll) implies 



where Q is the quintessence field, in general a linear combination of the moduli. The 
typical age of the universe being 10^*^ years, ()0.2|) suggests a very small variation of 
Q and hence of the physical constants (including the cosmological constant), well 
within the observational limits since big bang nucleosynthesis (BBN). This presents 
an alternative scenario where most of the evolution of the scalar fields occurs before 
BBN, while since then all the scalar fields are effectively frozen We discuss in 
details some of the consequences and observational bounds of having "rolling" moduli 
fields. 

However this is not the only possible scenario. Depending on the various parame- 
ters our potential can have one or two minima. The first one, which is always present, 
corresponds to the symmetric unsquashed state where the internal manifold preserves 
all the G isometries, while the second one can arise due the presence of the extra fiux 
and its depth depends on the value of the fiux. Thus it is possible that the depth 

more difficult question to answer. One should mention that in brane world scenarios, specially in 
the six dimensional version jSj, there is reason to believe that the quantum corrections may be of 
the right size. 
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of this minimum is so large that once our universe gets stuck it remains there. This 
leads to a different scenario where all the moduli are fixed! Of course in this case 
the fine tuning issues of the classical or bare cosmological constant remains unsolved. 
What could be really interesting in these situations is that even though the internal 
manifold is now frozen it could have been dynamic before, with possible application 
to inflation ^1 . Imagine the universe starting out at the "shallow" symmetric 
minimum and then making a transition, either to the second minimum or towards 
the "runaway behaviour". As explained in [1^1 this would correspond to a gauge 
symmetry breaking mechanism, G — > if, in the effective four dimensional world^; 
typically the mass of the broken gauge bosons corresponds to the scale of Grand Uni- 
fied theories (GUTs) fHlll]- During this phase transition, if the inflationary slow roll 
conditions [THj are satisfied around the top of the potential barrier, then our universe 
will experience an accelerated growth or inflation. Let us now briefly focus on some 
of the important details of the analysis and results. 

We developed general and physically intuitive techniques to study the dynamics 
governed by a sum of exponential potentials where the exponents are arbitrary linear 
combinations of two fields''. These techniques are based on a rich body of literature 
[THl US 1201 exponential potentials and in priciple could be extended to more 
than two fields. Also, recently a more formal analysis of critical points have been 
performed The analysis relies on finding the asymptotic behaviour of the potential as 
one moves towards infinity along any direction on the "moduli-plane" . Since we have 
only exponentials, ultimately only a single exponent dominates; the potential either 
grows exponentially or falls to zero. In this way one can map the entire moduli-plane 
into regions where each exponent dominates (see figEI)- Since the dynamics of a 
single exponential is exactly solvable one can construct qualitative trajectories of 
the "particle" as it passes through different regions dominated by different exponents. 
What is remarkable is that we find there are two ultimate late time attractor states of 
the particle (one when the shape, and the other when a linear combination of the shape 
and the radion is stabilized) and both of these solutions correspond to an effective 
exponent which is less than -\/2 signalling a phase of acceleration or quintessence. On 
the other hand all other exponents are larger than a/2 corresponding to radiation- 
matter domination^! 

Let us try to summarize the dynamics. We investigate three possible initial con- 
ditions depending upon which cosmology may vary: (I-l) the universe starts at the 
symmetric minimum, (1-2) the universe starts away from minimum with large (^ Mp) 
potential energies and as it rolls down it is able to completely avoid the symmet- 
ric minimum, (1-3) although the universe starts away from the minimum, there are 
mechanisms which depending upon the parameters can suck the particle close to the 
minimum after which the dynamics is similar to that of (I-l). We find that there are 

^In the braneworld interpretation squashing of the internal manifold will not be related to any 
gauge symmetry breaking in the brane because the bulk gauge fields are different from the ones 
living on the brane. 

^See |17| for a recent and more formal analysis of critical points in such systems. 
®For radiation domination one really needs the exponent to be larger than and this is also 
acheived for several values of the various dimensions in our model. 
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three possible late time states: (L-1) quintessence phase when the shape is stable, 
(L-2), another quintessence phase where a combination of the shape and the volume 
is frozen, and (L-3) the universe is stabilized at the second minimum. The interme- 
diate dynamics is typically non-accelerating (radiation-matter dominated) except a 
transitional phase where the particle comes out of the symmetric minimum in (I-l) 
and (1-3), and this can be accompanied by inflation. The aim of this paper is not to 
come up with "The Scenario" where everything works but to rather motivate how dif- 
ferent phases of cosmology: inflation, radiation-matter domination and quintessence 
can arise from moduli dynamics of SUGRA/String theory. Also, the paper develops 
analytical techniques to analyze complicated potentials involving exponentials that is 
expected in unified gravitational theories. These should be useful in further research 
to evolve a realistic cosmological history. Finally, the reader is warned that some of 
the analysis may appear too technical and detailed but considering the generic nature 
of the potential for SUGRA models the authors feel that it is important to study the 
dynamics in details. 

The paper is organized as follows: In section 1, we introduce our higher dimen- 
sional SUGRA type model with a brief review of pertinent group theoretical concepts. 
In section 2, we perform a consistent dimensional reduction including the fluxes. In 
section 3, the effective action and couplings to radiation is derived. In section 4, tech- 
niques to analyse dynamics with exponential potentials are developed and apphed 
to understand the general features of the cosmological evolution of the moduli. In 
section 5, we try to take steps towards realistic cosmological scenarios with empha- 
sis on inflation and quintessence. In particular, we study the various observational 
constraints on inflationary and quintessence cosmology. Finally, we summarize our 
results and indicate the challenges that need to be overcome. 



As in [iniE] we consider our universe to be a semi-direct product, M^+i (S> G, where 
Md_|_i is the D + 1-dimensional observational universe and G, a Lie group manifold^, 
serves as the Kaluza-Klein internal space [21] • Let us first review the Lie group 
geometry in brief. 

Geometry of Lie groups: A Lie group element g can be parameterised as 



where To G Q, the Lie algebra corresponding to the Lie group G and x"(2/™) ^-^e 



some given functions of the coordinates charting the Lie group manifold. The Lie 

^Although for simplicity we choose to work with a group manifold, most of the analysis depends 
only on the various dimensions entering into the problem and hence should hold at least qualitatively 
for more general homogeneous spaces like the coset spaces. 
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group generators To satisfy the usual commutation relations: 

[To,To] = Coo^To (1.4) 

O 

where Coo^ are the structure constants of the Lie group. With each of the generators 

ab 

To, one can associate a left and a right invariant vector field Co and Co respectively. 

We are interested in metrics which are at least invariant under the action of the 
right invariant vector fields. The components of such metrics are constants if we 
choose the left invariant vector fields Co, defined by 



mo m / a\— i 

Co = Co ao Co = I e o I 

a a in' a ^ m ' 

and 

g-^dog = eo^To , (1.5) 

as the local veilbien basis. A special case of the left invariant metric is the bivariant 
which is invariant under both the left and the right invariant vectors. The Killing 
metric given by 

g^, = -Coo^Coo^ (1.6) 

is an example of such a metric. Further, the Killing metric satisfies Einstein's field 
equations 

o 

-Roo =A goo (1.7) 

ab ab ^ ' 

and hence is consistent with its usual identification as Kaluza-Klein vacuum, the 

o 

constant A being referred to as the internal curvature. Contrary to this picture of 
an internal manifold frozen in its maximally symmetric Killing metric, we treat it 
as dynamic. In particular, it is interesting to study whether the manifold makes a 
transition from the Gl^Gr Killing metric to a ® -ff_R-invariant "squashed" metric, 
thereby effecting a gauge symmetry breaking from Gr Hr in four dimensions |l5j. 
The metric in this case looks like 




Here and from now on we refer group quantities by a circle (o) while that of the 
Cbset space G/H and the Subgroup H with (c) and (s) respectively. Sometimes we 
may omit the symbols when it is self-evident. We will also assume the group G to be 
semi-simple and the coset decomposition to be reductive 

CJ = (1.9) 

ab 

In ()1.8|) \I'^ corresponds to the "size" of the internal manifold and 0^ is the "squash- 
ing" parameter. Note, that we are in a broken phase when 7^ 1 and thus, as noted 
in [16t .4^. symmetry breaking can take place via two kinds of transition: (a) can 



5 



make a transition from a symmetric vacuum = 1 to a non-symmetric vacuum 
at = Bq 7^ 1 indicating a Higg's like mechanism, and (b) B^ can keep evolving 
very slowly towards infinity much like a quintessence field, effecting what we call a 
quintessential transition! 

Action and the Ansatz: As explained in the introduction, our aim is to try to 
connect the cosmology and particle physics that we observe in our four dimensional 
world with unified higher dimensional theories, like Supergravity/M-theory. Typi- 
cally the bosonic sector of such theories consists of the metric, some form fields and 
occassionally the dilaton. Thus, our starting point is the higher dimensional action 
given by 

where (f) is the dilaton field and the field strengths -F„j^'s are nj forms, / = 1 . . .M 
and we have also included a cosmological constant term. 

To obtain an effective four- dimensional theory one has to perform a consistent 
dimensional reduction jS]. One then ends up with Einstein's theory of gravity coupled 
to several "reduced" form fields and a bunch of scalars corresponding to the moduli of 
the internal manifold. Since one can usually interpret the dilaton as the radial moduli 
corresponding to a one dimensional circular compactification, we won't include it 
seperately in our discussion from here on. Also, although both in the action and in 
the field strengths one can have Chern-Simons-like terms, they can be ignored for the 
ansatz that we make. 

Since we want to study the symmetry breaking mechanism we have to include the 
squashing mode, while the size or the breathing mode has to be included for con- 
sistency. Including the shape or the squashing field has the additional advantage of 
being able to produce non-trivial curvature potentials^" which can enrich the cosmol- 
ogy that one obtains. In the Kaluza-Klein reduction scheme we thus have to upgrade 
\1/ ^ \E'(x) and B B(x). We will denote the coordinates charting the observable 

universe M^i+i by while x"^ will be used to collectively denote {x™,?/™}. As is 
perhaps evident by now, we are using "hatted", (^), quantities to refer to objects 
corresponding to the full higher dimensional manifold. 

Although an expression of the metric of the form (jl.8|) is physically clarifying, 
technically it is more convenient to include the scalars in the vielbein. We choose to 
parameterise the group element as 

9 = exp{ xHynT^Jexp{xhz^)T.J (1.11) 

^°In a cosmological context such potentials were first studied for specific coset spaces in j221 
focussing on cosmological instanton and domain wall solutions. In ^ squashing in arbitrary compact 
group manifolds was discussed while in cosmological solutions, where the internal manifold is a 
product space including having a hyperbolic geometry, has been studied. More recently in [24, three 
dimensional group manifolds, both compact and noncompact, have been studied in details. 
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The ansatz for the full higher dimensional vielbein is then given by 



and 



V 








\ 







;i.i2) 



:i.l3) 



The flat-metric is then just a constant 



9: 



ab 



Vab 

g^o 



:i.i4) 



We did not include the vectors in the ansatz ()1.12lll.l^ because we are only interested 
in the vacuum dynamics and the vectors appear as fluctuations around the vacuum 
metric. 

For the fluxes, several different choices of ansatz are possible depending upon the 
degree of the form field. Note that we are also not interested in the dynamics (or 
fluctuations) of the form fields but rather the condensate or vacuum expectation value 
which typically depends on the moduli fields once one solves the field equations that 
follow from [22] 



1 



{n-l)\ 



2•••P/^ 







and the Bianchi identity 



dF = 



:i.i5) 



:i.i6) 



Thus the F"^ term in the action gives us an additional potential term for scalars. 

Now, depending upon the exact internal manifold and the non-trivial forms that 
it can support one can have non-zero fluxes in different directions. However, in our 
model some of the choices appear naturally: 
1. Flux along G: 



2. Flux along H: 



mi. ..mo mi. ..mo 

D D 



mi...ms mi...m, 



Equipped with the ansatz we can now proceed to perform a consistent dimensional 
reduction. 
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2 CONSISTENT DIMENSIONAL REDUCTION 

Effective Action from Gravity: First let us look at the Einstein-Hilbert term in 
the higher dimensional action ()1.10|) . To compute the scalar curvature corresponding 
to the ansatz ()1.12lll.l'i|) it is convenient to use differential forms. One starts with 
the basis 1-forms 

= dx%; (2.1) 

For (jl.l2j) the 1-forms are given by 

cD^ = <ifQj (2.2) 

One then proceeds to compute the Ricci tensor via the spin connections and the 
curvature 2-forms in the following manner: The spin connections u"'^ are uniquely 
defined via 

rfb" + cD\AcD^= (2.3) 
The curvature 2-forms can be computed from the spin connections: 

n^^=duf^ + u;\Au;% (2.4) 

The coefficients of the Riemannian tensor can now be read off from the curvature 
2-forms 

^% = R\s:p'^^' (2-5) 

Without going into the details here, |1| we enumerate the Ricci tensor obtained from 
the usual contraction of the Remannian tensor 

Rm = R\i (2-6) 
After some simplifications one obtains 

Rab = Rab- b ^-V5(e„^)- b 0-^6(6,6)- b ^ Q'^ C ^a"^ e,)Q 

R,. = ^,c[-{$-in^ + {b ~l)^-\d^)^+ b '^-^e-^da-^d^'Q} 

ab ab 

2 \ 2 2 

= ^.j-i^-iQ^ + e-^ne + (b -1)^-2(9^)2 + (b -i)e-\dQ? 

ab ab 

+ {b + f))^"^e-i9a^9"0} + ^^~^Q\l -k) + Q-^k}] (2.7) 
Here we have assumed^^ that 

k gll^ CJ-^CJ = kCsJ^Cso^ (2.8) 

"■b ab ac hd ac m 



These assumptions are not true for an arbitrary subgroup if of G. 
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and 



ac bd ac hd 



(2.9) 



Note k and /c' are not independent group theoretical quantities but are related by 

(2.10) 



1-k' = '^{l-k) 



c 

D 



For a symmetric coset decomposition, i.e. Ccc'^ = we have 



k' = 0^k = l- 



D 

s 

2D 



(2.11) 



Another interesting case is when H = U{1) or a product of f/(l)'s. Then 



k = 0^k' = 1 



2D 

c 

D 



(2.12) 



We are ready to compute the scalar curvature that we need to compute the 
Einstein-Hilbert term in the action (|l.l(Jj) : 



S 



D 



Now 
or, 

R = R- 



R = g^'R^^ = g'^'Rab + g'^'R.j + g'^'R 



ab 



at 



2D — + 2D — + D {D -l)^-^+ D {D 







^2 



2 , 2D{D+l) 



(2.13) 
(2.14) 



4^2 



c s 



D +2 (1 - A;)- D{1- k)Q' + k D jr^ 



(2.15) 



Since R is independent of the group coordinates one can perform the integration 
over the group in the action ()2.13|) which essentially just yields a volume factor 



Vg = j d^y e 

Thus we have our effective D + 1-dimensional action 
Vg 



(2.16) 



5. 



grav 



OS 1 



16nG 

where we have identified 



e-'^^Q^R 



IQnG 



o s A/f r OS 



e-^^^e^i? = ^ / e-^^^e^i? (2.17) 



M2 



2 167rG^+i IGTrGg 
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f2.18) 



Mp being the four dimensional reduced Planck mass, 10^^ Gev. For most part of 
the analysis we work in natural units where Mp = 1. It is useful to perform some 
integration by parts. The simplified action looks like 



(2.19) 



where we have defined the Kinetic and Potential like terms^^ for the scalar fields as 



K 



D {D -1)^^+ D {D -1)^0^ + 2DiD-V 



(2.20) 



and 



V 



4^2 



iD +2D{1- k))+ D{1- k)Q^ ~ ^ ^ 



(2.21) 



One can also compute the effective action term corresponding to the higher di- 
mensional cosmological constant 



Sr. 



-2A 



dx 



(2.22) 



The Fluxes: Next, let us solve the Bianchi identity ()l.l(i|l and the field equations 
fll.l5|) for the field strengths. For simplicity we assume that only one of the n-forms 
discussed above is turned on. 

1. Flux along G |2|: The Bianchi identity essentially tells us 

F? , =c{y)eo o (2.23) 

mi. ..mo mi...m„ ^ ' 

D 



D 



Then 



ooo OO lo 00 

D,mi...mo o , , mi.. .mo i o-J-o2 , . o , , mi. ..mo /r^r^,\ 

F D = (g^)^c{y)e D = — -g e (y) c {y)e o (2.24) 

■q,2DQ2D 



D 



Plugging ^TM in (frTH|l we find 

c(y)=c.e(y)-i (2.25) 



with the flux energy reading 



f [D + l)\f' (2.26) 



^^Note that the coefScients in front of the kinetic terms only depend on the various dimensions 
and is therefore completely general for any internal manifold. The coefficients of the potential terms 
however will vary as one passes from group manifolds to more general manifolds, although the 
functional dependence on O and ^' remains the same. 
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Sf = --MII Id^^'xe-'^^ (2.27) 



The effective action is then given by 

4 J -qjDQD 

o — 1 o 

where we have absorbed 9 in the definition of c. The above flux contribution is well 
known in literature but we now come to a more non-trivial flux contribution along 
the subgroup H which as we shall see later, is consistent to turn on under certain 
conditions. 

2. Flux along H: As usual the Bianchi identity implies 

F? , =c(z)es o (2.28) 

mi...ms mi...ms ^ ' 

D D 



Then 



D,mi...ms s , , mi...ms gl s-^g , , m^,, ^ 



F D = {g c {z)e D = - e (z) g c {z)e 



mi...mu+i . 



Now 

^ = e-\x) {y)\flD{z) e'' {z)^^e^ 

c 

where D{z) is the determinant of the adjoint representation of H in G/H, Dc^{h{z)) 

and e (y) is the determinant of the vielbein along the coset space G/H. This is 
possible because in our choice of co-ordinate system the vielbein decomposes 

so that 

e^hy,z) =eJ{y)Df{hiz)) 

As we shall soon discover, for consistency one requires D{z) to be a constant, d. 
Plugging the metric and the field strength now in the field equations (|1.15|) one finds 
that it is satisfied if 

c (z) =c . (z) (2.29) 

The vacuum flux energy reads 

„ s2 s s-l 1 

F' =c D\ 9 — r (2.30) 

and the effective action is then given by 

1 9 r \lj{D-2D) 

- - c / rf^+^x e-'^-^ (2.31) 



4^ 



s s-l 

where again we have redefined c to absorb 9 . 



QD 



Consistency of the Truncation: In the previous subsection we obtained the di- 
mensionally reduced field theoretic action for our model. It is important that we 
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check the consistency of our ansatz i.e. the solutions that one obtains by varying 
the effective action (ll.lUj) should indeed be solutions of the full higher dimensional 
Einstein's equations. In ||4j it was shown that the truncation is consistent when we 
only have the Einstein-Hilbert term. More specifically it was found that 



\-qfDQD-l 50 



s 



^^^) g^. = b Gs. (2.32) 



and 

i - .\ ^^r ) % = -M|e- b (2.33) 

Consistency of the ansatz (|1.12lll.l'ij) is now obvious^^ from (|2.32j) and (|2.33p . 

We now want to investigate whether the consistency is violated when we introduce 
the fluxes. To do that we need the field equations for the higher dimensional metric 
corresponding to the full higher dimensional action p.lO|) : 

- - y ^(n/F-r p f/^-^"/ - -g^F^ ) + A^^^ = (2.34) 

ab 2 ^-^ 77,/! ac2...c„j b '2^ab rijJ ^ ab ^ ' 

Now, as in the case with no flux, there is no field redefinition involved in the ordinary 
space-time component of the metric and hence we only need to worry about the field 
equations along the internal dimensions. Specializing to a single flux we thus have 

Goo - -^(nFo^ ^ ^^...cn _ \g^,Fl) + A^oo = Goo - - = (2.35) 

ab 2n\ ac2...c„ b 2 ab ^' ab ab ab ab ^ ' 

In order for the consistency arguement to hold we need to show^^ 
and 



gs. = Mle-^ D Tss (2.36) 



5'S'fiux 1 5'S'fiux \ ^ ,^2-1 



g,, = Mle-^ b Tec (2.37) 

ah P ah ^ ' 



Let US verify this case by case: 
1. Flux along G: It is easy to see 



1 1 o ^C2...Co 1 o2o-l / 1 \ 

-— DFo Fo o = -cg — ^ g.o 

2 JJ\ ac2...c, b 2 \^,2DQ2D '^^ 



Then from ()2.35p we have 



1 o2o-l / 

-eg - 



ab 4 \^,2DQ2D J '^^ 



^■^Since g™" = g""\ i.e. there has been no field redefinition involving the four dimensional part 
of the metric it is clear that field equations corresponding to the variation of the four dimensional 
metric is automatically consistent. 

^''It is easy to check that the cosmological term satisfies the analogous equations. 
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Now let us compute the variation of the reduced action: 

SSfiux 1 ,^2 -1 Ao2o-l 1 



A« DC9 „ . (2.38) 



and 



Then 



Similarly, 



-M',e-'Dcg ^ (2.39) 



\'. 'tn" ) h = >p<=" ^-^'^^ = *^p<="' ^ 



- ^ h ^ (D -D)cg = My D T . 



2. Flux along if: As in the earher cases it is straight- forward to compute Too. We 
find 

1 520-1 / 1 \ 

T,,^- 'eg ^ _ g^, (2.40) 

ab 4 \^,2DQ2D ) 

and 

1 s2o-l / 1 \ 

T.. = - c 5 ^ g.s (2.41) 

4 \^2D02D / 

The variations are given by 



o s 



- Mle-^ DC 9 — ^ (2.42) 



50 4 ^ QD+i 
and 



5* 



--M'e-\D -2D) eg (2.43) 



Then again we have 

y-qtDQD-i 6k3 J ""-b 4 ^ ^,2bQ2b ^ "■b 

Similarly, 



)-l dO "-b 4 ^ ,Tf2D«2D ^ " 



hl = -AV -'^D + D)cg 



\^D-10£) 5^ ^D0£)-1 5<d ) <^b 4 ^ ^2L)02£) 

= M^'e-^ b Tec 

P ab 

Thus we have shown that we can have a consistent dimensional reduction even with 
flux turned on along G or H. Also, since in general the degree of the flux-form is 
different in the two different cases, the field equations and consequently the consis- 
tency does not interfere with each other, and one can imagine all the fluxes turned 
on simultaneously. 
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3 EFFECTIVE ACTIONS AND COUPLINGS 



In the previous section we performed a consistent dimensional reduction of a higher 
dimensional supergravity type action containing the pure Einstein-Hilbert term along 
with a cosmological constant and fluxes. To analyze the cosmological and particle 
physics aspect of the model, it is however convenient to perform conformal rescalings 
so that the Einstein-Hilbert term is canonical in four dimensions. 

The Conformal Transformations: Let us first redefine the moduli fields to bring 
the kinetic terms ()2.20j) in the usual form: 



^ = e^; and 9 = e'' 
In the new variables the gravitational action looks like 



gravj 



where 



K 



D {D -l){di)f+ D {D -l){def + 2D{D ~l)daijd^e 



and 



'grav ^ 

One can also rewrite the cosmological and the flux terms. 

M? 



and 



S. 



flux 



1 

4^ 



o2 

c e 



Di^-De_^ ^{D-2D),p-De 



We now perform the conformal rescalings: 



D , D 



The total effective action then reads 

Seff = iSgj-av + '5'cos + 'S'fiux ~ 2 j '^^^^'^ ^ ^[R + K ~ V], 

where 

1 



K 



D-2 



D{D- l){d^Y+ D{D+D- l){dey + 2D{D- l)dai^&'e 



(3.44) 

(3.45) 

(3.46) 
(3.47) 

(3.48) 
(3.49) 

(3.50) 
(3.51) 

(3.52) 



2 
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2 



with 



s 

D 



(3.53) 



o 

D 



P 



7 



and 



D-2 

o 

D.D 
D-2 

s c 

D D+D 
D-2 

b {D+ b) 



(3.54) 
D-2 

Diagonalization and Normalization: In the earher subsection we obtained the 
complete effective action for the moduli scalars coupled to gravity. For convenience 
in the cosmological analysis, it is useful to diagonalize the kinetic terms in the action. 
Such a diagonalization can be acheived by choosing the volume {v) and the shape {6) 
as variables 

V — %l) -\- v6 

After diagonalization we have 



1 



D-2 



b{D- i){dvf + ^{b +D - i)(de)^ 

D 



(3.55) 



while the potential in terms of this new variable looks like 
1 



-2v 



o2 

(3.56) 
(3.57) 



where 



5' — 5 — 



c s 

DD 

o 

D 



(3.58) 



For cosmological analysis it is also convenient to rescale the fields 



V — > KyV and 9 — > KqO 



(3.59) 
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with 

b{D-i) bb 

K^ = \ and Ke = (3.60) 

\ D-2 ^ D 

so that the coefficient in front of the kinetic terms is normahzed to the canonical 
value of 1/2. This essentially has the effect of rescaling the exponents for the fields v 
and 9 by and K^^ respectively. 

Couplings: We have successfully obtained the classical effective potential (scalar 
interactions) for our moduli fields. However they will in general also couple to matter 
and radiation. These couplings will depend on whether the matter-radiation fields 
are residing on the brane [22] or on the bulk |2j. We focus on the radiation^^ term 
and obtain the couplings for the two different cases: 

1. Bulk Radiation: Fluctuations of the off-diagonal part of the full higher dimensional 
metric can be identified as four dimensional gauge fields. These fields will in general 
couple to the moduli fields and the coupling can be obtained following the usual 
dimensional reduction procedure One finds that after the conformal rescalings 
the action for the gauge fields to be given by 

Sbui,,„ = -\j d^^'x e-^e^(^+^)F^ = / d^^'xe-'e'^^-'^'-^^'^ (3.61) 

s s 

where is the usual field strength for the unbroken (non-)abelian gauge fields A'^. 

2. Brane Radiation: In the brane- world picture the action for radiation localized on 
the brane is given by 

Sbrane,gauge = J d'^X ^-^(brane)F„„Fp,r''^"« (3-62) 

where T is the brane tension in Planckian units. Now ^(brane) is related to (^(brane) 
through conformal rescalings 



-^(brane) = ■^/-^((brane) A"(^+^) 
while 

and thus the A dependence completely cancels in the brane action for D = 4, which 
is only a consequence of the conformal invariance of Yang-Mill's action in four dimen- 
sions. Then ^ 

Sbrane, gauge '^'^ j ^ ^ (3.63) 

i.e. there is no direct coupling of the moduli fields to brane radiation term. 

Obtaining fcrmionic couplings are more complicated as the fermionic representations can be 
different depending upon the number of extra dimensions. Also, four dimensional fermions can 
come from higher dimensional gravitinos. We leave this topic for further study. 
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4 COSMOLOGICAL EVOLUTION 

Exponential Potentials and the Mapping Technique: Having obtained the 
effective potential let us try and analyze what it looks like in the v — 9 plane. To 
have a qualitative understanding of the potential we develop a "mapping technique" 
which works in general for a potential with sums of exponentials of two fields. We 
start by looking at how the potential behaves as one goes asymptotically to infinity 
in any direction, say 

9 = mv 

Since the potential is a sum of exponentials it is clear that as one proceeds to infinity 
along a particular direction only one of them, the one with the largest (or smallest) 
"effective exponent" will dominate as v ^ +oo (or — oo). Thus in the v — 9 plane one 
can map out regions where the different exponents dominate. Moreover, depending 
upon the sign of the dominant effective exponent and the coefficient in front of the 
dominating exponential potential, the potential will either increase ("wall"), decrease 
("fall") or asymptotically approach zero ("plane"). The map therefore provides a 
qualitative picture of the dynamics. Since there are six exponentials in our model one 
has to compare six different effective exponents for a given slope, a.s v ±oo: 

Xi = -2 + 2(1 + u)m 

Xa = -2 + 2z/m 
X3 = -2-2(1 - u)m 
X4 = -2a 
X5 = -2(3 

c s 

Xe = -2(7 + m^) (4.1) 
D 

Also note 

P > 'J > a while 1 > z/ > (4.2) 

Comparing the different exponents in 1)4.11) one can obtain ranges of m where a par- 
ticular exponent dominates. For example, for m > and f ^ 00 it is clear that 
Xi > X2 > X3. While 

1 — a 

X4 > Xi ^ -2a > -2 + 2(1 + z/)m ^ m < 

i.e. X4 overtakes Xi when m becomes less than (1 — a)/(l + u). Proceeding in this 
manner one can indeed map out the entire v — 9 plane. Computing the sign of the 
dominant exponent also tells us how the potential looks as one goes towards infinity. 
Consider the Xi exponent: 

1 

Xi > ^ -2 + 2(1 + iy)m > ^ m > 

^ ^ 1 + z/ 
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Thus when m > 1/(1 + z/) the potential increases exponentially (the coefficient in 
front being positive) so that we hit a wall. However, when m < 1/(1 + z/) then 
the exponent is negative, so that the potential asymtotically falls off to zero and 
the same thing happens as Xi gives way to X4. Proceeding in this manner once we 
know the entire v — 6 map, it is easy to understand the dynamics, specially when one 
is sufficiently away from the origin^^. Then essentially the "particle" rolls along a 
single (the dominant) exponential potential and such behaviour has been extensively 
studied [TH] in the literature. 

If the dominating exponential term is given by 

^ Ve'""^'^^ (4.3) 

with V being the coefficient in front then the particle moves along a constant direction 
given by 

u = -^^{av + p9) (4.4) 
under the influence of an effective exponent 

"eff 



v/«2 + f32 (45) 



Now it is known ^Hj that a^ff < indicates a scalar field dominated regime of 
accelerated expansion while agff > corresponds to matter-radiation domination 
where the scalar field tracks matter-radiation densities. So let us find out whether 
the exponents we have support acceleration or not. To obtain bounds on the effective 
exponents, we of course rescale the exponents using ()3.(i()|l . Playing around with the 
dependence of the various dimensions in the exponents we find^^ 

Xieff > V2 



Xaeff > V2 

Xseff > V2 
V2 > X4eff > 

/e 



3V^> X4eff >3\^>V2 
V 5 

Xeeff > 2 (4.6) 

It is clear that the only region which generically supports acceleration is when X4 
dominates. This is not surprising given its origin in the higher dimensional cosmolog- 
ical constant but what is important in our model is that, we do not need to start with 



-'^^The origin gets its special status because we assume that all the mass parameters are of order 1, 
i.e all the potentials are roughly of equal strength near the origin. If this is not the case, the origin 
essentially gets shifted, but the analysis still holds. 

^^There may be stronger bounds but these are sufficient for our purpose. 
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a "small" cosmo logical constant. One should be cautioned that this behaviour is only 
expected at late times if and when the particle is well "inside" the region when one 
can approximate the potential by a single exponent. The particle can have interesting 
transient behaviour when it is moving from a region dominated by one exponent 
to a region dominated by a different exponent. In our model, one example would 
be when the particle enters X4-region from Xi-region. Under certain conditions the 
particle may also follow a different attractor solution which is a compromise between 
the two attractor solutions corresponding to the two competing adjacent exponents. 
Such situations will arise in our model with interesting consequences and we will 
elaborate on them later. To begin with let us try and understand the generic picture 
with the help of the v — 9 map. 

Quintessence or Moduli Stabilization: To understand the effect of the fluxes we 
first look at the potential when we have not turned on the fiux along the subgroup. 
In this case (see fig. P) the minimum in the potential corresponds to the symmetric 
internal manifold (gauge symmetry unbroken) which makes it interesting for studying 
possible symmetry breaking mechanisms jU . 

We note that the moduli in the lower half plane is not stabilized and in particular 
the potential falls to —oo in the region bounded by the straight lines 

6 = V and a = v 

1-u l-v 

essentially because dominates in this region. Thus to perform meaningful cosmol- 
ogy one has to be able to avoid the "fall" which may be possible depending upon 
the initial conditions. Before we investigate this any further though, let us see what 
happens if we now turn on the fiux along the subgroup. 

The potential looks (see fig. Ej) similar to the previous case except that we find the 
shape moduli field 9 now has a stable potential, perhaps not surprisingly. Previously 
it has been argued in the context of string theory and Calabi-Yau manifolds how 
fiux should stabilize the shape moduli fields pp. Our model furnishes an explicit 
realization of the same, albeit in the context of an internal group manifold. 

Let us now focus on the cosmological evolution of the moduli fields. We will 
consider the most general case when all the coefficients of the exponents are non-zero. 
The evolution when some of the coefficients are zero (like when the subgroup fiux is 
not turned on) can mostly be derived from the general case, although interesting 
differences sometimes emerge and we will comment on them as we go along. One 
can get a good idea of the evolution by studying the v — 9 map (see fig. EI) but 
the evolution clearly depends on the initial conditions. The Initial conditions can 
be broadly classified under two categories: either (I-l) our universe starts from the 
symmetric minimum (see figure) and then may evolve to a broken phase, or (1-2) it 
starts out far away from the minimum in any particular direction. Among 1-2 cases 
for realistic cosmological scenarios, we consider the particle to be initially rolling down 
a wall {Xq,X^ or Xi). Let us consider these situations in a little more detail. 
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Figure 1: The potential with c = 0. 

I-l. There are two possibihties depending on how deep the hole (potential minimum) 
is. 

(a) The tunnelling probability may be so small that the particle remains trapped for 
the entire evolution of our universe thus far; there will be no interesting cosmology to 
be gained but at least we would have solved the moduli stabilization problem^^. Note 
a subtle but interesting difference when the extra flux is turned on: the minimum 
is no longer at the symmetric = point. It is slightly shifted indicating a small 
squashing. In other words if the hole is deep then the universe will be stuck in a 
slightly squashed or broken phase. As observed in H] the mass of the broken 
gauge bosons depend very sharply on the squashing and in order to produce a large 
hierarchy the squashing has to be very small, which may be achieved in this scenario. 

(b) If the hole is shallow then the particle may be able to get out of it by tunneling. 

""^^ Observe that the internal manifold will be trapped in its symmetric phase which often preserves 
some amount of supersymmetry 2 and hence has a chance to address the cosmological constant 
problem. 
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Figure 2: The potential with c 7^ 0. 



It is however well known ^3] that the collision of bubbles formed as the universe 
undergoes a phase transition through tunnelling produces inhomogeneities that are 
too large as compared to observation and hence we do not consider this possibility 
in any detail here. However, if the hole is really shallow so that a small velocity-kick 
originating from classical or quantum fluctuations may induce the phase transition 
then a realistic inflationary scenario may emerge as was considered in j^j. We should 
point out that when one introduces the additional flux, then for certain parameter 
ranges, the potential can develop a second minimum. If this second minimum is stable 
then the inflaton can oscillate around it to reheat the universe ^2] and eventually 
settle down offering us a scenario where at Late times all the moduli are stabilized 
(L-3). Alternatively, one can imagine a hybrid type inflationary model j^H] where as 
the inflaton approaches the minimum, it veers off in another "steeper" direction (and 
thereby ending inflation) and ultimately to either L-1 or L-2. This is reminiscent 
of the quintessential inflation models j2Z] where nevertheless the role of inflaton and 
quintessence field is played by different fields. We will study these in some detail in 
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c s 

D D 



Figure 3: The v — 9 map for A; > 0. The sohd radial hnes demarcate regions where different 
exponentials dominate, the arrows indicating the direction in which the particle moves. The 
dashed line corresponds to Qfine = constant. 

the next section. 

1-2. From the v — 9 map and the arrows in the map (see fig. EI) one sees that, in 
general, along the upper half plane one moves from the Xs-region onto the Xi-region 
(or starts from Xi region) and then finally to the X4-region. Note, since neither 
nor Xi supports acceleration, the universe while passing through these regions would 
evolve as in the radiation-matter dominated phase. This will eventually be followed 
by an accelerated phase in the X4-region which can continue for eternity and we will 
call this as the L-1 Late time behaviour. Because of the presence of the extra flux, 
the evolution in the lower half is no longer pathological. One typically moves from the 
Xs-region to the Xg-region (or starts from Xg itself) and then finally to the X4-region 
where we can have quintessence. Again, since Xq also does not support acceleration, 
the generic scenario is radiation-matter domination followed by quintessence. 

Thus, for 1-2 we seem to have a history of our universe consistent with observa- 
tion; radiation-matter domination followed by a phase of quintessential acceleration 
which presumably we are observing today. It should be mentioned that in [22] it was 
also observed that one could have late time acceleration phases for compact inter- 
nal manifold when relevant fluxes are turned on. However in order to have a viable 
quintessence scenario one has to also satisfy the cosmological observational bounds 
on the quintessence equation of state and the cosmological constant energy density. 
Further, if the quintessence field couples to matter/radiation, observational bounds 
coming from time variation of various coupling constants and fifth force experiments 
jTU] need to be satisfied and we will investigate all these in the next section. 

With the help of the maps we have so far discussed qualitatively some likely evo- 
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lution of the universe. One thing we should mention is that it is not necessary to 
assume that our current universe hes in the X4-region, although this is the only region 
which supports acceleration. The reason is transient acceleration may potentially be 
able to explain the observed phenomenon of quintessence [21]. For example, we can 
be caught up in the transition between the Xs-region and the Xg-region, or between 
say, the Xs-region and the Xi-region. One can also have other viable late time accel- 
erating phases as we will see in the next subsection. 

Dynamics with Two Exponential Potentials: We have seen that for a wide 
range of initial conditions SUGRA type models give rise to a radiation-matter dom- 
inated phase followed by quintessence. However one can also obtain an early phase 
of acceleration (inflation) from the dynamics near the well. As mentioned earlier, if 
our universe starts from the symmetric minimum and the potential well is shallow 
then the universe could inflate |4 while getting out of the well via an initial velocity 
kick. Here we will find that it is not necessary for the universe to start from the well 
but that there are mechanisms available which can lead the particle towards the well 
even if one starts with initial condition 1-2. This can happen if the particle is trapped 
in a "trough" created by two competing exponents as we will illustrate. Thus, the 
I-2-type initial conditions further divides into two categories: either the particle fol- 
lows the original trajectory described earlier and we will continue to refer such initial 
conditions by 1-2, or the particle can be attracted towards the well, type 1-3 initial 
conditions. In this case its further dynamics proceeds in a manner similar to type I-l 
initial conditions, so in some sense 1-3 is a mix of I-l and 1-2. To understand how this 
happens we need to analyze the dynamics of the particle when there are two adjacent 
competing exponential potentials. As a bonus we will also find a second late time 
accelerating attractor solution! 

[ZQ E3 EZZ] 

i ii iii 

Figure 4: Three different late-time behaviours depending on the orientations in competing 
regions (see text for details). 

For simplicity, let us consider a potential containing two exponentials: 

2 

V = \4e"i^+^i^ + ^/se^^^+z^^e ^ ^ y,^x, . ^ (4.7) 

i=l 

As our earlier analysis showed, the v — 9 map gets divided into two regions, where the 
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two different exponents dominate (see ffgE)). Now, depending upon the direction of 
motion in these two regions which can be computed from (j4.4|) . three different kinds 
of dynamics or attractor behaviour is possible as illustrated by fig. El In case (i) 
the arrows point towards each other so that a "trough" is created between the two 
regions. So no matter where the particle initially starts from it will ultimately end 
up in the trough and continue to move exactly along the boundary demarkating the 
two competing exponential regions 

Mcomb = ±(/?2 - Pl)v T ("2 " ^1)6 (4.8) 

The particle evolves as if governed by an effective exponent given by 

"1/^2 - "2/?! 



X, 



eff 



y'(a2-«i)2 + (/32-A)^ 



(4.9) 



Thus it may happen that even though the individual exponents are too large to 
support inflation the combined exponent ()4.9|) may, and this was already noted in 
some special cases in [T^ . 

Coming back to the different cases, in (ii) we see that the arrows are pointing 
in opposite directions. Such a situation implies an unstable ridge between the two 
regions. So there are two late time attractor solutions. The particle either slides to 
the Xi-region or X2-region depending on the initial conditions and evolves as if the 
other exponent {X2 and Xi respectively) is absent. 

In (iii) one of the arrows points towards the other region but the second arrow 
points towards its own region. The dynamics is again simple. The second arrow 
is the late time attractor. The particle may evolve in the first region for a while 
according to the exponent Xi but eventually it encounters X2-region, after which it 
makes a transition to the X2 late time attractor solution. The dynamics during such 
a transition may also be quite interesting as pointed out in 



Phases of Acceleration : Let us investigate when in our model one can obtain 
troughs of type (i). In order to do this first we need to figure out the directions (uj) 
of the trajectories in the different regions and then check when we get type (i) for the 
different adjacent regions that are relevant. Here is the list: 

Ml = 2{v-{i + u)e) 

U2 = -2{v- vO) 

M3 = -2{v + (1 - v)e) 

U4 = 2av 
U5 = 2i3v 

ue = 2{-fv + 6'e) (4.10) 

Clearly, the "velocities" ui and uq imply a transition of type (iii) between Xi — X4 
and Xq — X4 (see fig. E}, and thus the late time behaviour is dominated by X4. Let 
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us now compare the velocities Ui and Uq with the slopes demarcating the X5 — Xi 
and X5 — Xg region respectively. One finds that we have type (iii) situation if the 
various dimensions satisfy 



2 



< 2 (4.11) 



o s s 2 o 

{D- D) D (D +2] 

o 1 s o 2 

(D +1 D) D 
and 

+ < 2 (4.12) 

(b-i)b 

respectively and the dynamics proceed as in 1-2. Otherwise, depending upon the 
initial conditions as the particle nears the X5 — Xi or X5 — Xg boundary it will be 
sucked into the X5 — Xi^^ or X5 — Xg combined track (type i) and the particle will 
enter the "near- well" region; we refer to this case as 1-3. 

We further observe that near the well, the exponents X2 and X3 kick in, which 
can be instrumental in creating an attractor behaviour in conjunction with other 
exponents. Suppose k > 0. Then as one approaches inward towards the well in the 
Xg-region, the X3 exponent become strong as can be seen from 1)4.11) . A trough or 
a valley is created, sandwiched between the Xg and so to speak, the Xs-region, as is 
also evident in the picture (fig. Ej) of the potential. One can compute the effective 
exponent (Xag) in this valley. 



X. 



36 



<v^ (4.13) 



\ b+2 



Thus this can indeed support acceleration. This phase of acceleration gives us an 
alternative late time quintessence phase, L-2. We observe that if the higher dimen- 
sional cosmological constant is absent, then this trough will actually be the generic 
late time attractor whereas if k = then L-2 is absent. 

To summarize, we have three initial conditions 1-1,2,3 and two late time quintessence 
evolution, L-1,2, plus a stable L-3 late time configuration. For 1-2 L-1,2 the dy- 
namics is clear; radiation domination followed by quintessence. However for 1-1,3 — *• 
L-1,2, 3 there is an intermediate evolution in the vicinity of the well about which we 
about which we try to shed some light in the next section. 



5 REALISTIC SCENARIOS AND CONSTRAINTS 

Interesting Cases: We now have a qualitative understanding of the different tra- 
jectories that the particle may follow, as well as its cosmological significance, but the 
possibilities are numerous. So here let us focus on the likely and interesting ones. 

o 

^^We note in passing that for large number of extra dimensions, £)> 30, the X^ — Xi track supports 
acceleration and can therefore account for a phase of inflation. 
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First of all, let us assume that the well is shallow enough so that the particle can 
roll out of it. Second, we realize that depending on the initial conditions and the 
various parameters of the potential, the particle/inflaton, after crossing the potential 
barrier, can enter a region dominated either by Xi, X^/Xq or X4. In the first two 
cases inflation can end as neither exponents support acceleration; however, it doesn't 
if the infiaton directly goes over to the X4-region. Hence the realistic scenarios can 
be summarized as: 



Case 



I.e. 



Evolution 



Phases 



1. 



I-l 



Xi X4 



E^R^L-1 



2a. 
2b. 
2c. 



I-l 
I-l 
I-l 



W- 
W- 



W2 
X3 



36 



W^X3,X6^X4 



E^R^L-2 
E^R^L-1 



3a. 
3b. 
3c. 



1-2 
1-2 
1-2 



X5 
X5 



Xi X4 

Xq X4 
Xq X36 



R-^R-^L-1 
R^R^L-1 
R^R^L-2 



4a. 
4b. 



1-3 
1-3 



X5 
X5 



Xi5 ^1. or 2. 
Xes ^1. or 2. 



R^R^l. or 2. 
R^R^l. or 2. 



where W2 labels a stable second minimum, E denotes a possible early phase of inflation 
and R corresponds to radiation-matter domination. As is evident from the table, most 
of these scenarios possess a late time acceleration phase which may be able to explain 
the observations today, while some of them admit an early phase of acceleration or 
inflation. In the next subsections, we investigate in details whether these phases of 
acceleration can indeed be compatible with the various observational bounds coming 
from early and late time cosmology as well as particle physics. 

Constraints from Inflation: For an inflationary scenario to be successful it should 
be able to produce at least 50-60 efoldings to explain the standard cosmological issues 
presented by the flatness and horizon problems ^H]- The inflation scale should be 
such that it generates CMB fluctuations with the right amplitude (5h) ^H] 

5h ~ 10"^ (5.14) 

Now, to answer these questions comprehensively one has to perform a detailed 
numerical analysis involving evolution of both fields, but one can get significant insight 
into the physics by employing approximate analytical methods as we will show. We 
start by "freezing" a linear combination of the fields j3], i.e. assume that the infiaton 
moves along a constant slope in the v — 6 plane: 

cr = f + m^"'^^^ = constant (5.15) 

and for analytic simplicity we choose 

"^infl^O (5.16) 
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The subsequent evolution of the inflaton can also be approximated by straight line 
trajectories the slope of which is determined by the specific regions it is traversing 
(refer to table). A simplified picture of evolution would thus look like fig. El Of course 
in reality, these straight line evolutions will be joined by smooth curves. 



m=-(p-l)/(l+v) 




m=l/(l+ V) 



m=(l- a)/(l+v) 




Figure 5: The v — map for /c = 0. The dashed Hnes correspond to straight line approxi- 
mations to the more realistic (solid Une) evolution. 

We start our analysis with the potential ()3.56p which can be rewritten as 



-2d 



D 



o2 

+2Ae-2"^ - -b + -e'^/^" 
4 2 

where we have also made the approximation that near the well, the Xq contribution 
is negligible. It is known that the minimum of the potential V{9) occurs at 

^min = and V{e = 0) = 

so that the total potential is minimized by minimizing V!nin(^)- One can compute this 
approximately. In general for a potential of the form 



with 



a\> a2> OLz 



(5.17) 
(5.18) 
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/ V2a2 \ «i-a2 

and 

vJ^-i] flips'; ^ 



With some algebraic manipulations we find in our case 



D+2 



Thus if Vfiux is large compared to typical Kurv, the typical hierarchy required being 
10^, then Kin is much smaller as compared to Mp and produces the right size of CMB 

fiuctuations ()5.14|1 . For example say D= 2 and Vfiux ~ 10"^. Then 



4 ^= 5 



Next let us try to estimate the number of e-foldings. For this we have to compute 
the slow roll parameters 

q2 

eH = ?> -<1 (5.21) 

and 

riH = « 1 (5.22) 

rlu 

First observe that 

V{d) = Kin + \e-'''{V{e) + K(^n.in, 'c)e-^''') = \e-''- {2 + V{d) + V^e'^''') (5.23) 

where we have reintroduced the Xg potential as it plays an important role away from 
^ = 0. V{9) looks like a Mexican hat. Now infiation can occur as 9 rolls down from 
the maxima to the second minimum along which the potential can be approximated 
by a cubic potential, the derivative of which has to vanish at the maximum {—9x) as 
well as the minimum {—9m)'- 

V'{9) = c{9 + 9,){9 + 9m) (5.24) 

In terms of 

C = -(^ + ^m) (5.25) 
V{0 = Kn + c{\c' + ^AC^) (5.26) 

where 

Kn = V{-9m) and /\ = 9m-9^ (5.27) 
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With ()5.26|) one can compute the total number of e-foldings as 



^M+ATa+A/'g (5.28) 

—Q ^ —A is where the inflation starts while — Ce ~ is where the inflation ends. 
Notice As is small and the main contributions to the e-foldings come from A/i and 
A/2. If we do not assume any fine tuning of the initial conditions the logarithms are 
pa2 so that 

A-«4l| + ^ (5.29) 

Now all we have to do is to identify Vm, c, 6m and Ox from the actual potential ()5.23|) . 
Note that Ve is the last to become active, so that 6x is still given to a good approxi- 
mation by V{6) and one knows this value exactly [1^1 E] 

s c 

e''^ = ^-4^ (5.30) 
2D + D 

Next, we realize that the second minimum is produced essentially because of the 
competition between X3 and Xq so that using and ()5.2H|1 we get 

o 

Om = ^^ ln7^; n=^ (5.31) 

2DiD -1) ^6 

and 

Vm = (^b +2 - k{b -l)7^'/(^"'^) (5.32) 

71 is basically like a measure of the ratio between Kurv and Vfi^^^n- Further one can 
compute V" exactly at this minimum 

D 

which gives us c. Using ()5.30H5.3^ then, one obtains an expression for A/" in terms 
of the parameter TZ, a plot for which looks like fig. IHl First, observe that even the 
minimum e-folding that one gets is around 10 which is already of the right order 
that is required to explain the cosmological problems [Jj. Second, note that one can 
get large number of e-foldings in two regions of small and large TZ although caution 
should be made in the small TZ region as after a point the second minimum is above 
the first symmetric minimum and the analysis will break down. 

The case of small TZ inflation when most of the e-foldings come from TVi signals 
a hybrid inflationary model [2E]- Note here Vm 7^ and therefore inflation does not 
end through 6 oscillation but rather as the particle approaches the minimum it starts 
to roll off in another direction. In the context of our model a natural choice would be 
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Figure 6: No. of e-foldings vs. TZ 



the volume v direction or perhaps a second squashing mode which is responsible for 
later quintessence! For large TZ inflation, most of the e-foldings come from A/2 with 
the inflaton eventually settling down at the minimum (L-3). In this case it is clear 
that we cannot solve the cosmological constant problem. TZ has to be fine tuned to 
make Vm = 0. In general one has to add a four dimensional cosmological constant by 
hand to make it vanish. 



Constraints from Quintessence: From the table we see that there are really two 
possible late time attractor solutions, the X4-region and the X^q trough. In either 
case the particle moves along a constant slope given by 



iriq = = for X4-region 



and 



rrin 



1 — a 
l-u 



for Xse-trough 



(5.34) 



(5.35) 



()5.35p follows by applying ()4.8p . During this evolution, the combination of fields 

= 9 — TJigV (5.36) 

remains "frozen" and the particle rolls along the orthogonal direction Q effectively 
along an exponential potential 



(5.37) 



where the "exponent" a„ is given by 



2a 



2D 

'ib+2 



for X4-region 



(5.38) 
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and 

c 



0!l2 = "36 



2 D 

for Xae-trough (5.39) 



\ b+2 



It is now clear how a modest rolling of the quintessence field can produce the famed 
120 orders of mismatch between Planck mass and the cosmological energy density [7] 

^ 10-120 ^ y^^^ ^ (5 40) 

p 



From inspection of ()5.38|5.39|1 . it is clear that aq is smallest when D= 1 ^ ctg ~ ^2/3 
while it is always less than so that 



V2> aq> ^2/3 ~ 0.8 (5.41) 

To generate the hierarchy ()5.40|) then, we just need a rolling of Q from say Qo to Qc 
given by 

Ag=|Qc-Qo|~ — 

150 > AQ > 85 (5.42) 

Let us look at the bounds on equation of state parameter uj coming from obser- 
vational cosmology [Zj: 

u < -0.78 (5.43) 

Now, for exponential potentials the exact late time attractor solutions are known 
For ttq < \/2 this corresponds to a power-law phase of acceleration 



2 



a{t) ~ t'^^ (5.44) 
with an equation of state parameter Uq given by 

LJq = -l + (5.45) 

Then, compatibility with ()5.43|) seems to tell us 

aq < 0.8 (5.46) 

which is just possible when D= 1 ()5.45|) . Unfortunately for compact group manifolds 

C 

D= 1 is not possible although for coset space reductions^" and non-compact internal 

c 

manifolds (in the brane-world scenario 28j) it is possible to have D= 1. We find that 
t)=2^ aq = l while £), £)= 3 ^ ~ 1.1 



^°An interesting example will be squashing a 2-spliere in a 6-dimensional SUGRA context, which 
has become specially popular as a brane world model (Si- 
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These clearly does not satisfy the current bound ()5.43|) . Fortunately this is not the end 
of the story. Recently it was discovered that although the late time attractor solution 
for an exponential potential satisfies (|5.45p . before reaching the attractor solution 
there is a transient phase of acceleration where the equation of state varies |24J. In 
particular it was shown in that this transient phase of acceleration can explain 
the current phase of acceleration with uj close to -1, even though ujq is not. Taking 
this point of view let us now investigate the further constraints coming from violation 
of equivalence principle and time variation of fine structure constant experiments. 

Imposing ()5.36p one obtains an effective Lagrangian which now depends on a single 
scalar field Q: 

~l~ -^scalar -^rad 



iSeff — d X 



with 



and 



(5.47) 



:scalar = -^(5g)2-A,e-'« (5.48) 

/:rad = --^e-'^^F' (5.49) 



The constant parameters Of and Xg depend on Q±. We have used (j5.36p and normal- 
ized Q so that the kinetic terms are canonical 



Q = Kv; K=^K^ + mlK^ (5.50) 

We now investigate the bounds on the coupling exponent of Q to the electro-magnetic 
F"^ term. Note that in our model for radiation in the bulk, the fine structure constant 
is time dependent and is given by 

Of = aoe'"^ (5.51) 

Then 

I — I = |a|Q < 10-V~' (5.52) 

Of 

where the bound comes from atomic clock experiments For solutions ()5.44|) it 
is known that the kinetic energy is a fraction of the scalar potential energy which is 
identified with the effective cosmological constant. This gives us a handle to estimate 
the Q variation 

K. E. = (D P. E. ; (13 = ^ 

l-UJq 

Recent CMB data (23 seem to suggest Uq ~ —1.02 ±0.1 using which we find 

^ g ^ 10"^^ - lO^^V"^ (5.53) 

This implies a bound on a 

\(t\ < 10"^ - 10"^ (5.54) 
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Let us compute a for L-1 and L-2 late time solutions. For bulk radiation from ()3.61|) 
one gets in general 

. = ^1^ + '^ (5.55) 
K 

In particular when = — 1/(1 — z/), cr = and the bound is obviously satisfied. 
From fig. El it is clear that the L-2 slope is closer to the o" = line as compared to 
the L-1 slope and hence the situation with the bounds 15.551 is improved. However it 
is not close enough. Substituting mLi and in ()5.55|) we find 

(Tli = K^> \f2 and ^ < crL2 = < ^2 (5.56) 

respectively. Thus unless the uOq parameter is significantly closer to -1, ctli cannot 
satisfy ()5.54|) while (Tl2's case is marginal. For brane radiation no such problem arises 
as it is not coupled to the moduli fields. This finding however is consistent with the 
earlier realization that perhaps our current universe is in a transient quintessence 
phase rather than having already reached the late time attractor solution. 

Finally we come to the most crucial observational bound concerning variation of 
Of on cosmic scales. Assuming to be small one finds 

^ ^ -aAQ ~ -aQAt (5.57) 

Of 

where At is say the time elapsed since BBN. For BBN it is known that ^ < IQ-^. 
On the other hand we have seen that typically to be consistent with quintessence 
cosmology |A(5| ~ 100. Thus the first part of the equality gives a very strong bound 
on cr < 10^3 - 10^"^. Not only do our late time cr's violate this bound but also it is 
hard to see how transient quintessence can explain such small cr's without fine-tuning. 
However, this has two possible resolutions. Either (a) one has to invent a mechanism 
such that the quintessence trajectory has a late time attractor slope very close to 
(or possibly equal to) —1/(1 — i') when a = 0. This may be possible to achieve by 
introducing other flux terms or considering potentials coming from brane gas [n] etc. 
Or (b) if |AQ| is small! How can this be? Well, all we need is for Aobbn to be small. 
If we look at the second part of the equation ()5.57|) . then substituting the typical age 
of the universe since BBN, At ~ 10^°|/r, we find |AQ| ~ 10^^! Then even for ex ~ 1, 
is well within the observational bounds. Of course this is only an approximate 

argument and Q need not be this small through out but (b) does seem to offer a viable 
alternative. What this means of course is that most of the hierarchy in the potential 
has to come from early (pre-big bang) evolution of Q during inflation, reheating etc. 
Crudely, the evolution of the particle can be represented as in fig. El where it starts 
out with an initial pre-BBN slope but then slowly changes (or is still changing) to 
the late time quintessence slope. 

It may now seem that we have solved the problem posed by ()5.57|) . but in the 
bulk case not quite! The reason is that in our potential ()3.56|) X3 is identical to the 
radiation coupling exponent ()3.61|) modulo a — sign. For the current value of Of this 
means that X3 provides a large contribution to the cosmological constant and neither 
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resolutions (a,b) can work without resorting to fine tuning of tlie parameters in order 
to cancel the X3 contributions, but this is precisely what quintessence models aim 
to avoid! The only way to avoid these fine tunings is to ensure k = when the 
potential term is absent. Remarkably, one finds that k, which is determined by 
group theory, vanishes whenever if is a U{1) or a product of ?7(l)'s. This is nice if 
one eventually wants to make a connection to symmetry breaking in Standard Model. 

Finally, one can look at the various null experiments on fifth force ^21- In general 
the quintessence field will also couple to matter 

S'n.atter = / d'^x V^xif+ ime^'^)x (5.58) 

and this will lead to violation of equivalence principle JU] ■ However typically bounds 
obtained from these experiments come in the combination of /i^ and /icx: 

/icr < 10"^° (5.59) 

and this only tells us that for the values of a in our model 

< 10"^ (5.60) 

and we leave the study of fermionic couplings to future research. 

6 SUMMARY AND FUTURE RESEARCH 

To summarize, we have derived potentials for the volume and the squashing mode in a 
typical SUGRA model. The potential gets contributions from the internal curvature, 
the flux, and the higher dimensional cosmo logical constant. It turns out generically 
to be a sum of six exponentials with rich features which we explored to understand 
the corresponding cosmological phases. We find that at late times the moduli could 
either be stabilized around a minimum (L-3), or one of them (a linear combination 
of shape and size) could be frozen while the other evolves towards ibinfinity (L- 
1,2). In this latter case one obtains quintessence as the external universe undergoes 
accelerated expansion for these late time attractor solutions. Prior to reaching these 
late time solutions, the universe in general undergoes radiation-matter domination as 
the moduli potential is steep (dominating effective exponents being larger than a/2). 
Thus, for a class of initial conditions our potential predicts the history of our universe 
to be radiation-matter — > quintessence (acceleration). On top of these phases our 
model can also accommodate an early phase of acceleration or inflation if and when 
the internal manifold makes a transition from the symmetric (unsquashed) minimum 
towards the asymmetric (squashed) minimum. In the original Kaluza-Klein picture 
where the four dimensional gauge fields originate from the higher dimensional (bulk) 
graviton, this phenomena also corresponds to a gauge symmetry breaking mechanism. 
For brane gauge fields no such interpretation is possible. 

To understand how these various phases of cosmology emerge, both qualitatively 
and quantitatively, we developed techniques to determine the moduli evolution under 
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the influence of exponential potentials. These techniques work for any number of 
exponentials with two fields and perhaps can be generalized to more than two fields. 
Hence, one could add more fiuxes or introduce brane-gas etc, and study their effect on 
the evolution/stabilization. Finally, we should note that although we specialized to 
group internal manifolds, most of the results (including the exponents) only depend 
on the various dimensions in the problem and hence are completely general to compact 
internal manifolds. Some specific details may change because the coefficient in front 
of the exponentials can vary and it would be interesting to look at coset spaces, 
specially those whose isometry group is similar to the gauge group of the Standard 
model or GUTs. 

The analysis and the results so far are encouraging, but several problems still 
remain before we can envisage realistic particle physics and cosmology. First let us 
focus on the quintessence aspect. From our study of several constraints coming from 
variation of fine structure constant and the fifth force, several fine tuning issues persist 
with regards to the radiation coupling exponent o. This is enhanced if one assumes 
universal coupling or Brans-Dicke type theory. There are essentially three ways to 
avoid these constraints: (i) one invents a mechanism so that the linear combination 
that is stabilized corresponds precisely to a = 0. Although we came close, we didn't 
really succeed. There is still scope of adding other fiuxes and/or introducing brane 
gas etc. which can do the trick but again this will only work for certain specific 
dimensions! (ii) one finds a mechanism like chameleon mechanism to ameliorate 
bounds on a . (iii) show that most of the scalar evolution that occurs precedes BBN, 
while since (or sometime before it) BBN the scalars have really slowed down so that 
they are effectively frozen. 

Next let us look at the infiationary scenario. The main problem associated with 
the infiationary scenario in our model is to have a graceful exit. If the second mini- 
mum is really deep then infiation can end in the usual way with oscillations reheat- 
ing the universe. However, in this case one cannot address the issue of the small 
cosmological constant as one has to fine tune the parameters to ensure that the po- 
tential energy at the minimum is zero (or close to it). The other choice is to have 
a quintessential-inflation type scenario where the quintessence fleld and the inflaton 
are nevertheless separate flelds. As the inflaton approaches the minimum it veers off 
in another (steeper) direction as in hybrid inflationary models except that this new 
direction (or fleld) plays the dual role of a quintessence fleld later on! This is a nice 
picture particularly because our model naturally contains two flelds (for example Q 
could be the inflaton, while v plays the dual role) except that it doesn't quite work. 
First, notice that quintessence has a chance to work only if A; = 0, on the other hand 
we found a second minimum (and flattish potentials) only when A; > 0. Second, if we 
assume /c > for successful inflation, the particle veers on to the trough which 
also supports acceleration. Similarly, if we assume = the particle veers off along 
V direction under the influence of X4 exponent, again signalling an accelerated phase. 
Hence ending inflation in these scenarios is a nontrivial task. We believe one can ad- 
dress this issue in two different ways, (i) As pointed out in [21] the particle does not 
immediately fall into the late time attractor solutions (X36 or X4) but rather there 
is a transitional phase of radiation-matter domination. One can investigate whether 
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such a phase can account for the radiation era in our universe, (ii) A more attractive 
scenario would be to invoke double squashing! Two successive transitions, one with 
k > supporting inflation, followed by a = transition taking care of the graceful 
exit problem along with later quintessence. 

This second scenario becomes even more appealing in the light of our common 
understanding of particle physics. Our universe is supposed to have undergone two 
gauge symmetry breakings - first, a symmetry breaking of a GUT gauge group to the 
Standard Model gauge group and second the Electro-weak symmetry breaking. Note 
for the second transition A; = 0, a is involved which is precisely what we need for 
the quintessence to work! The challenge of course now is to come up with the right 
internal manifold. If one wants to accomodate all the symmetries within the seven 
extra dimensions one has to look into coset spaces whose consistent truncations are 
slightly more difficult to study. On the other hand, even as a toy exercise, if one wishes 
to first develop a group manifold model, one has to generalize our constructions. This 
is because the Standard Model gauge group is not simple, not even semi-simple. This 
would involve including more squashing fields corresponding to the different simple 
components in order to have a consistent truncation and we leave this exercise for 
future research. 

To conclude several challenges still remain but perhaps we are just beginning to 
understand the implications of extra dimensional dynamics. 
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